Subwavelength imaging with opaque left-handed nonlinear lens 
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We introduce the concept of subwavelength imaging with an opaque nonlinear left-handed lens 
by generating the second-harmonic field. We consider a slab of composite left-handed metamaterial 
with quadratic nonlinear response and show that such a flat lens can form, under certain conditions, 
an image of the second-harmonic field of the source being opaque at the fundamental frequency. 
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One of the most unique properties of the recently 
demonstrated left-handed metamaterials, i.e. materials 
with simultaneously negative real parts of dielectric per- 
mittivity and magnetic permeability pj , is their ability to 
focus electromagnetic waves by a flat slab of the material, 
the property which makes these materials quite different 
from the conventional optical lenses with the positive re- 
fractive index needed to have curved surfaces to form an 
image. Recently, Pendry argued that a slab of a loss- 
less left-handed material with e = fi = — 1 should behave 
like a perfect lens enabling to obtain an ideal image of a 
point source through the amplification of the evanescent 
components of the field. 

While recent experimental demonstrations confirmed 
the main features of negative refraction of the left-handed 
materials 0, , the question of near-perfect imaging by 
flat lens and near-field focusing still remain highly con- 
troversial [5j , and it severely constrained because of large 
dissipation and dispersion of metamaterials. Neverthe- 
less, numerical studies indicate [|J that nearly-perfect 
imaging should be expected even under realistic condi- 
tions when both dispersion and losses of the left-handed 
composites are taken into account. More importantly, 
some of the properties of the left-handed materials, such 
as negative refraction, have been predicted and re- 
cently demonstrated experimentally [8| in photonic crys- 
tals, which are inhomogeneous periodic structures with 
a lattice constant comparable to the wavelength. This 
negative refraction allows considerable control over light 
propagation, and it opens up doors for new applications 
of microstructured materials and the extension of the ba- 
sic concepts from microwaves to optical frequencies. 

Until now, all properties of flat lenses and subwave- 
length imaging have been studied only for linear waves. 
However, it has been already noticed that the left-handed 
metamaterials may possess quite complicated nonlinear 
magnetic response [9j, or properties of such materials 
can be altered by inserting diodes in the resonators of 
the composite structure thus making the response of the 
entire material nonlinear |10|. In this Letter, we make 
one step forward into this direction and study the effects 
of the second-harmonic generation and subwavelength 
imaging by a nonlinear left-handed lens. In particular, we 



analyze the imaging properties of a slab of metamaterial 
with quadratic nonlinear response and demonstrate, both 
analytically and numerically, that such a slab can form 
an image of the second-harmonic field of the source being 
opaque at the fundamental frequency. This can happen 
under certain conditions, which include Pendry's condi- 
tions of a perfect lens satisfied for the second-harmonic 
field, e(2w) = /i(2w) = —1, leading to the conditions for 
e(ui) and fj,(ui) being of the opposite signs at the funda- 
mental frequency to. More importantly, for the case of 
two sources we show that the resolution of such a nonlin- 
ear left-handed flat lens can be made indeed better than 
the radiation wavelength. 

We consider a lens in the form of a layer of left-handed 
metamaterial with the thickness D, as shown schemat- 
ically in Fig. 2] We assume that the metamaterial is a 
three-dimensional composite structure made of wires and 
split-ring resonators (SRRs) in the form of a cubic lat- 
tice. When the lattice period d is much smaller than 
the radiation wavelength A, i.e. d -C A, this composite 
structure can be described within the effective-medium 
approximation, and it can be characterized by dielectric 
permittivity and magnetic permeability which, for the 
specific structure and in the linear regime can be derived 
consistently and can be written in the form 



c(w) = 1 - 



/Lt(cj) = 1 + 
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where 



= yj 2irc 2 /d 2 In (d/r) is the effective plasma 

frequency, w = ^oyi^—F), is the eigenfrequency 
of an isolated SRR, F is the volume density of SRRs, 
7 e and 7 m are the damping coefficients, to is the fre- 
quency of the external electromagnetic field, r is the 
wire radius, and c is the speed of light. In the frequency 
range where the real parts of e and fi are both negative 
and for 7 e ,7m <C ui, such a composite structure demon- 
strates left-handed transmission, whereas for lo < ujq, it 
is opaque because the signs of e and \i are opposite. 
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FIG. 1: Schematic of the problem. Electromagnetic waves 
emitted by a source (located at z — —z s ) are reflected from 
an opaque slab of a left-handed material. Inside the slab, the 
exponentially decaying field at the fundamental frequency ui 
generates the second- harmonic field at 2u>, which penetrates 
through the slab creating an image at the second harmonics. 



In order to achieve a nonlinear response in such a struc- 
ture, we assume that each SRR includes a nonlinear el- 
ement, e.g. a diode inserted in the SRR [lfj. If the 
diode has an asymmetric current-voltage characteristics, 
the unit cell does not possess a center of symmetry with 
respect to the field direction, and the resulting nonlinear 
response should include, in particular, the second har- 
monic of the source frequency. 

Our idea is to satisfy the perfect-lens conditions for the 
second-harmonic waves and, therefore, we take )i(2uj) = 
e(2w) = —1. From Eqs. QJ, © we can find that these 
conditions can be satisfied provided 

e(w) = -7, n(w) = (3-F)/{3-2F). (3) 

For this choice of the material parameters, the slab is 
opaque at the fundamental frequency w, and the waves do 
not penetrate into it. However, the nonlinear quadratic 
response of the material can allow the process of the 
second-harmonic generation. Since the material is trans- 
parent at the frequency 2uj, we expect that the second- 
harmonic field can propagate through the slab creating 
an image of the source behind the flat lens. 

Using the so-called undepleted pump approximation, 
we can obtain the equation for the TM-polarized second- 
harmonic field Hy (x, z) inside the slab, which has the 
form well-known in the theory of the second-harmonic 
generation (see e.g. Ref. 

K^)H^=- 1 -^^M^\ (4) 



AH, 



» 



where A is the Laplacian acting in the space (x, z), 

K 2 (2oj) = 4fcge(2a;)/^(2u;), and is the nonlinear 

magnetization of the unit volume of the metamaterial at 
the frequency 2u>, which appears due to the nonlinear 
magnetic momentum of SRR, 
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where 



c 3 d 3 U c R d uj 2 



-1 -i- 



Hy (x, z) is the spatial distribution of the magnetic field 
at the fundamental frequency in the slab, a is the ra- 
dius of the resonator rings, Ra is the differential resis- 
tance of the diode at zero voltage, fco = w/c, and U c 
is the diode parameter defined from the current-voltage 
characteristics of the diode which we take in the form 
I = I (e u ' u * - 1). The right-hand side of Eq. © van- 
ishes outside the metamaterial slab. 

Applying the Fourier transform along the x direc- 
tion, we obtain the equation for the function Hy' = 

By (k x ,z) in terms of H y ^ — Hy(k x , z) and G = 
G(k x ,z) which are the Fourier transforms of H^' (x, z), 
H^^x, z) and [H^ (x, z)] 2 , respectively: 



d 2 H 



(2u0 



dz 2 



[K 2 (2u;) - 4fc 2 ] H^ = r)G, 



(6) 



where r\ = 127rfcge(2w)x(a>). Using the convolution theo- 
rem, we express the function G through the spectrum of 
the magnetic field at the fundamental frequency Hy^ in 
the form 



G = 



H^\k' x ,z)H^\k x -k' x ,z)dk' x . (7) 



Within the framework of the undepleted pump approxi- 
mation, the function H^\k x , z) can be found as a solu- 
tion of the linear problem describing the electromagnetic 
field at the fundamental frequency transmitted into the 
left-handed slab, 



2n\e 



D s (u>, k x ) 



Z ie koK2Z - Z 2 e 



k K 2 (2D-z) 



(8) 

where S(j) is the spectral function of the source at 
the fundamental frequency located at the distance z s 

= V 1 ~ 7 2 > K 2 = 
Ki ± «K2/e(w), and 



from the left-handed slab, 



Zl,2 



sji 2 - e(uj)[i(uj), 7 = k x /k , 

D(u, k x ) = Z 2 — Z 2 exp (2koK2D). The spectral function 
5(7) includes both fast propagating (7 < 1) and slow 
evanescent (7 > 1) spatial harmonics. Possible distor- 
tions of the second-harmonic image can be caused by the 
pre-exponential factor in Eq. ©, and the main effect is 
due to the pole singularity defined by the equation 



D(v,k x ) =0, 



(9) 



that characterizes the resonant excitation of surface po- 
laritons which are known to limit the resolution of a left- 
handed lens. In the case of a thick slab, i.e. k^^D 1, 
the resonant wavenumber of the surface waves can be 
found in the form 



7 2 

/ sp 



e 2 H - 1 



(10) 



Substituting the explicit expressions for e(u>) and n(u>) 
from Eq. |J3J| into Eq. (|10|l . we obtain a simple estimate for 
an expected resolution limit of the nonlinear left-handed 
lens in terms of the critical (limiting) wavenumber, 



7Um 



1 + 



7(3 -2F)' 



(11) 



However, the existence of this critical wavenumber does 
not necessarily limits the lens resolution and, in reality, 
the effect of surface waves on the imaging properties of 
the nonlinear lens depends on the efficiency of their ex- 
citation by each particular source. 

Analytical solution of the problem for the spatial spec- 
trum of the second-harmonic field transmitted through 
the left-handed slab can be obtained for narrow enough 
spectrum of the source, i.e., when the width of the source 
spectrum at the fundamental frequency does not exceed 
the value j c , where 



7c 2 « KuOmMI. 



(12) 



Then, we can use the impedance boundary conditions for 
the fundamental field at the interface between vacuum 
and the metamaterial slab at z — 0. Subsequent numeri- 
cal results indicate that this approximation remains valid 
provided j c > 7u m . 

To solve the problem analytically, we assume that the 
wave at the fundamental frequency u> penetrates inside 
the slab on a distance (the skin layer) much smaller 
tha n the slab thickness D, i.e. D5 1, where 5 = 
koy/—e(u))fj,(u>). Taking into account the actual values 
of e(u>) and /z(u>) ©, one can see that the fundamental 
frequency penetration depth, ~ (27rJ) , does not exceed 
A/ 17. Then, Eq. © can be re- written in the form 



dz 



" + [K 2 (2uj) - 4fc V] 



rye 



-28z 



A)(7), 



where ^0(7) = £(Y)£(7 ~ l')dl', and $(7) = 

e Sz H^ (j, z) does not depend on z. As a result, the 
general solution for the second harmonic inside the slab 
can be presented in the form 



H y 2u) {l, z) = C ie 2koK2Z + C 2 e- 2ka " 



-2Si 



, (13) 



where G\ t % are two constants which should be determined 
from the boundary conditions, and 



1 



^0(7) 



(14) 



We should satisfy the continuity of the tangential com- 
ponents of the magnetic H y 2u ^ and electric fields 
at the interfaces between air and the metamaterial slab, 
i.e. at z = and z — D. As a result, we obtain the sec- 
ond harmonic of the magnetic field behind the slab (for 
z > D) where the image is expected to form, 



H (2u 
V 



'(7 z) = --(l- l K2e ^ ) 
l7 ' ' 2\ W(2w) 



Q 2k K!(2D-z) 



(15) 
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FIG. 2: Numerical results for imaging by a nonlinear left- 
handed lens. Shown are the intensities of the fundamental 
field at the source location (solid) and the second-harmonic 
field at the image plane for different cases (stars, dashed and 
dotted lines), (a) Wide beam (the width is A) generated by a 
single source, D = A, z a = A/2, (b) Narrow beam (the width 
is A/10) generated by a single source; dashed-the image for 
D = A/10, z s = A/20; dotted - the image for D = 0.3 A, 
z s — 0.15 A. (c) Imaging by two sources separated by the 
distance A/5; dashed - the image for D = A/10 and z s = 
0.03 A; dotted - the image for D — 0.3 A and z s = A/5. 



For wide beams with narrow spectra, Eq. (|15fl can be 
rewritten in the form 



H^H-i z) = 



V 



e{2uj) 



uj) - e(u)(i(w)] 
-1 



1 - i 



exp{2k oKl (2D-z-z s )} / S(r/')S{j - j')dj'. (16) 



Thus, the squared field at the fundamental frequency 
acts as an effective source of the second-harmonic field 
and, as a result, the image of the squared field is re- 
produced by the nonlinear left-handed lens. This image 
appears at the point z lm = D — z s , and this result co- 
incides with the corresponding result for the linear lens 
discussed earlier 0. 
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FIG. 3: Intensity of the fundamental [(a,c)] and second- 
harmonic [(b,d)] beams (in units of the wavelength) for the 
problem of the second-harmonic generation and imaging by 
a nonlinear left-handed lens. (a,b) Wide beam (D = 5A, 
z s — 2.5A, and ao = A) and (c,d) narrow beam (D — A/5, 
z s = A/10, and ao = A/4). Solid lines mark the flat sur- 
faces of the nonlinear left-handed lens. Dashed lines show the 
predicted locations of the second-harmonic image. 



When the size of the source is comparable or less than 
the wavelength A of the fundamental-frequency wave, the 
problem cannot be treated analytically, and Eq. (0} has 
been solved numerically. In Figs. |2Ia-c), we present our 
numerical results for the intensity distribution of the in- 
cident beam at the source point and the field distribution 
of the second-harmonic beam at the image location, nor- 
malized to the field maxima. The actual amplitude of 
the electromagnetic field at the image location is lower 
than the amplitude of the source because of a finite effi- 
ciency of the process of the second-harmonic generation. 



For the objects with the spatial scale larger or equal to 
the radiation wavelength, the second-harmonic field pro- 
file coincides with the intensity of the fundamental field 
generated by the source, as shown in Fig.|3[a). 

However, when the source contains the spatial scales 
less than the radiation wavelength, the imaging proper- 
ties of the nonlinear lens depend strongly on the slab 
thickness D. As an example, in Fig. EJb) we show the 
results for the transmission of an incident Gaussian beam 
of the width A/10 which reproduces almost exactly the 
source profile at the image plane in the case of a thin 
lens (dashed line) but generates a strongly distorted im- 
age when the slab thickness becomes larger than a half 
of the wavelength A. Distortions appear as periodic vari- 
ation of the second-harmonic field being caused by ex- 
citation of surface waves. Intensity distribution of the 
magnetic field in the fundamental and second-harmonic 
fields are shown in Figs.|3Ja-d) for (a,b) large and (c,d) 
small (compared to the radiation wavelength) size of the 
source, respectively. 

Figure Etc) shows the numerical results for imaging 
of two sources that generate the Gaussian beams with 
the maxima separated by the distance A/5 in the trans- 
verse direction. Again, the image reproduces very well 
the source for a thinner lens, and therefore a thin nonlin- 
ear lens does provide a subwavelength resolution of the 
second-harmonic field. In contrast to the linear flat lens, 
the resolution of the nonlinear lens depends on the dis- 
tance z s between the source and the lens, and the optimal 
distance can be determined separately for each particular 
case. 

In conclusion, we have studied the second-harmonic 
generation and subwavelength imaging by a nonlin- 
ear left-handed lens, i.e. a slab of metamaterial with 
quadratic nonlinear response. We have demonstrated, 
both analytically and numerically, that such a slab can 
form an image of the second-harmonic field of the source 
being opaque at the fundamental frequency, with the res- 
olution that can be made indeed better than the radiation 
wavelength. 
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